Planck Spectroscopy and the Quantum Noise of Microwave Beam Splitters 
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We use a correlation function analysis of the field quadratures to characterize both the black 
body radiation emitted by a 50 load resistor and the quantum properties of two types of beam 
splitters in the microwave regime. To this end, we first study vacuum fluctuations as a function 
of frequency in a Planck spectroscopy experiment and then measure the covariance matrix of weak 
thermal states. Our results provide direct experimental evidence that vacuum fluctuations represent 
the fundamental minimum quantum noise added by a beam splitter to any given input signal. 
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At optical frequencies, single-photon detectors [1] and 
beam splitters are key ingredients for the successful de- 
velopment of atomic physics and quantum optics. These 
devices are crucial for the implementation of quantum ho- 
modyne tomography [2], quantum information process- 
ing and communication [3], as well as all-optical quan- 
tum computing [4]. The recent advent of circuit quan- 
tum electrodynamics (QED) [5-15] has paved the way for 
the generation of single photons in the microvifave (mw) 
regime [11, 15]. Despite the rapid advances in this pros- 
pering field, the availability of suitable mw photodetec- 
tors [16, 17] and well-characterized mw beam splitters is 
still at an early stage. However, we have recently shown 
that the use of low-noise cryogenic high electron mobil- 
ity transistor (HEMT) amplifiers represents a versatile 
and powerful approach for the analysis of both classi- 
cal and non-classical mw signals on a single photon level. 
Although the phase-insensitive HEMT amplifiers obscure 
the signal by adding random noise of typically 10-20 pho- 
tons at 5 GHz, they do not perturb the correlations of sig- 
nals opportunely split into two parts and then processed 
by two parallel amplification and detection chains. In 
such a setup, a correlation analysis allows for full state 
tomography of propagating quantum mw signals and the 
detector noise, simultaneously [18]. We note that HEMT 
amplifiers represent available "off-the-shelf" technology 
and offer flat gain over a broad frequency range of sev- 
eral GHz. Here, we present results of two experiments 
demonstrating the successful application of our setup to 
the characterization of weak thermal states. In a first 
experiment denoted as Planck spectroscopy we analyze 
the mw black body radiation emitted by a matched 50 f2 
load resistor as a function of temperature in the frequency 
regime 4.7 < uj/2tt < 7.1GHz. Besides confirming that 
the mean thermal photon number follows Bose-Einstein 
statistics [19, 20], our data directly show that the quan- 
tum crossover temperature Tcr shifts with frequency as 
= huj/2kB, providing an indirect measure of mw 



vacuum fluctuations with high fldelity. In a second ex- 
periment, we use weak thermal states for a detailed ex- 
perimental characterization of mw beam splitters at the 
quantum level. This task, which has not been accom- 
plished previously, is particularly important because mw 
beam splitters are key elements in a variety of quantum- 
optical experiments such as Mach-Zehnder and Hanbury 
Brown and Twiss interferometry [1, 19]. 

The experimental setup is shown in Fig. 1. As an ideal 
black-body source emitting thermal mw states [21], we 
use matched 50 fl loads whose temperature T can be var- 
ied between 20 and 350 mK and measured with a RuO 
thermometer. The temperature is measured by a RuO 
thermometer. The associated quantum voltage can be 
expressed as Vth = Vo{p^ +p), where Vq — 4BWi?o/iw/2, 
-Ro = 50 il, and p are bosonic creation and annihi- 
lation operators, and BW is the bandwidth. The ther- 
mal mw signal is fed to the input ports of a 3 dB mw 
beam splitter. We perform experiments on two different 
beam splitter realizations: a four-port 180° hybrid ring 
(HR) [cf. Fig. 1(b)] and a three-port Wilkinson power 
divider (WPD) [cf. Fig. 1(c)]. For the WPD, an inter- 
nal distributed resistor R^q shunting the output ports 
B and D provides isolation between those ports and 
impedance matching for port A. In addition, an exter- 
nal 50 load is attached to input port A. For the HR, 
matched 50 il loads are attached to both input ports 
A and C. The input /output relations of the HR are 
Vb = -?(VA + Vc)/v^and Vb = (-Va + Vc)/72 [22]. 
Remark that the WPD, although appearing to be a three- 
port device, has to be treated quantum mechanically as 
having an additional "hidden" internal fourth port C, 
see Fig. 1(c), ensuring energy conservation and commu- 
tation relations. In this case the input/output relations 
are Vb = (Va - Vc)/V2 and Vb = -i [Va + Vc)/V2. 
Regarding thermal noise, the internal resistor i?eq can 
be modeled as two equivalent matched 50 loads adding 
correlated thermal noise via the hidden port C only. If 
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FIG. 1: Schematics of the experimental setup (cf. main text), 
(a) Amphfication and detection chains, (b) 180° hybrid ring 
(HR). (c) Wilkinson power divider (WPD). For the WPD, 
port C represents a hidden internal port and 7?cq an internal 
distributed resistor, which can be modeled as two equivalent 
matched 50 Q, loads adding correlated thermal noise via the 
hidden port C only. 



the input signal at port A is large, this additional noise 
can be neglected and the WPD can be formally treated 
as a three-port device. 

As shown in Fig. 1(a), the output ports B and D of 
the beam splitter are connected to two symmetric am- 
plification and detection chains, each comprising the fol- 
lowing elements: (i) A cryogenic circulator at 500 mK 
with 21 dB isolation which prevents amplifier noise from 
leaking back to the mw source under study. Further- 
more, together with the isolation of two outputs of the 
beam splitters they avoid spurious correlations of the 
noise originating from the two amplification chains, (ii) 
A low- noise HEMT amplifier thermally anchored at 4.2 K 
with power gain Ghemt — 24 dB and noise temperature 
Themt ~ 6 ± 1 K. The noise added by the linear HEMT 
amplifiers is the dominating amplifier noise in each detec- 
tion channel and is expressed by x = + 0> where 
and ^ are bosonic creation and annihilation opera- 
tors [23]. (iii) A room-temperature amplifier and (iv) a 
mixer (M) and local oscillator (LO) to down convert the 
mw signal to the intermediate frequency (IF). The phase 
</jlo of one of the two LO signals obtained by means of 
a power divider (PD) can be varied by a phase shifter 
(PS), (v) An IF amplifier and (vi) a 7kHz to 26 MHz 
band pass filter. Using a double sideband receiver the 
total bandwidth BW in our experiment is twice the band- 
width of the bandpass filter. The down converted volt- 
age signal components si and S2 in the two chains are 



finally synchronously digitized by an acquisition board 
with 4 X 10*samples/s and 12 bit resolution. 

According to Nyquist [24], the black body radiation 
emitted by a resistor with resistance i?o within the band- 
width BW around the frequency lo = cjlo is given by 
(Vj.|)/i?o = 4BW(nth)^i^. Here, (nth) is the average 
thermal photon population. It has been shown that for 
conductors with a large number of electronic modes the 
statistics of the emitted photons is given by a Bose- 
Einstein distribution [21]. In this case the well known 
result {V^^)lRo = 4BW%^ coth(fia;/2fcBT) is obtained. 
The power emitted into a perfectly matched circuit with 
characteristic impedance Zq is reduced by a factor of 1 /4 
due to voltage division. Together with the beam splitter 
input /output relations of the HR, the signal components 
si and S2 are given by si — \/Gi{~ 



*/3^th + Xi) 

and S2 = y/G2{iaV^^ - ifiV^ + X2)- Here, d ~ 
G2 is the total power gain of the amplification chains, 
a — (3 = l/2-\/2, and xi X2 are the independent 
noise contributions of the amplifiers. Equivalently, for 
the WPD we obtain si = ^JUl{-iaV^^ + i^V^ + xi) 
and S2 = v/G^(-«aVfh - «/^Kh + X2)- We note that 



and xi,2 are classical realizations of the op- 



erators given above. By recording a large number of 
1/iS-long time traces (^ 10®), the auto- and cross- 
correlation functions i?ii(r) = (s*(t + t) Si{t)) /Zq = 
alsmc{BWT)/Zo and Rijir) = {s*{t + T)sj{t))/Za = 
sinc(BWT) cos(/jlo/-^o, respectively, can be calcu- 
lated — 1,2). Since {si{t)) — for thermal states, 
the auto- and cross-correlation functions are equal to the 
auto- variance Cii{T) = Ru^t) — (si)^ and cross- variance 
Cij{T) = Rijir) — {si){sj), respectively. Here, r is the 
time shift between two traces being correlated, af^ and 
afj are the variance and covariance of the voltage signals 
si and 32- 

We first discuss the Planck spectroscopy experiment. 
Here, we use only a single amplification chain and de- 
termine the auto-correlation function Rii{t) or R22{t). 
Exemplary, Fig. 2(a) shows the measured Rii{t) curve 
obtained for T = 30 mK using a WPD. A very simi- 
lar result is obtained for R22{t) [see Fig. 3(d)]. Fitting 
the data to Cii(r) allows us to extract the measurement 
bandwidth to BW ~ 52 MHz. Assuming that the signal 
contributions and V^^ due to the two load resistors 
and the noise Xi of the HEMT amplifier are independent, 
we simply can add up their variances and obtain Ru (0) = 
C,,(0) = al = {sj)/Z, ^ <g[(^a^+p^)lyl) + l^^j)]. 
With /iw/2fcBTHEMT <C 1 we can introduce a classical 
noise temperature TIjHEMT for the amplifiers and obtain 



Ru{Q) = G*BW 



"2" ^ 2k^ ^ 'i^BTi HEMT 



•(1) 



Here, G* — •jGi is the effective gain with Gi the total 
gain of the amplification chain and T'i*HEMT = 21,hemt/7 
is the amplifier effective noise temperature, representing 
the amplifier noise temperature relative to the input of 
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FIG. 2: Planck spectroscopy of thermal mw states at uj/2-k = 5.3 GHz using a WPD. (a) Auto-correlation function _Rii(r). The 
symbols represent the data, the line a fitted curve, (b) Temperature dependence of the variance i?ii(r=0) (Planck function). 
The symbols represent the data, the dotted and full lines are obtained by two and three parameter fits, respectively. The inset 
shows the data over a wider T range, (c) Planck spectroscopy. Contour plot of Rwij — 0) versus T for different Lj/2n. The 
data has been corrected for the frequency dependent amplifier gain, the offset due to vacuum and amplifier noise has been 
subtracted. The symbols and the dashed line mark the experimentally measured and theoretically expected quantum crossover 
temperatures, respectively. 



the WPD. For both, the HR and WPD, 7 = 4:{a'^ + P^) = 
1. 

Figure 2(b) shows the measured variance i?ii(0) as a 
function of T for a;/27r = 5.3 GHz in the case of a WPD. 
During the experiments, T was varied between approx- 
imately 20 and 350 mK by means of a resistive heater 
and continuously monitored with a RuO thermometer. 
As expected, the measured variance is close to a Planck 
function and nicely reproduces the quantum crossover at 
Tc: — haj/2kB- A fit of Eq. (1) to the data using 
and T]*HEMT ^ parameters yields G* ~ 90 dB and 
T^n^MT — 6K. Evidently, there are slight deviations be- 
tween the data and the two-parameter fit. This can be 
understood by taking into account that the effective elec- 
tronic temperature Toff of the load resistors at ports A 
and C may differ by a small amount 6T. Using ST as 
the third fitting parameter the solid line in Fig. 2(b) is 
obtained, demonstrating excellent agreement with the ex- 
perimental data. The ST values obtained by fitting the 
data are reasonably small and typically range between 
1 and 10 mK. The large bandwidth of the HEMT am- 
plifier allows us to perform equivalent measurements at 
several frequencies between 4.7 and 7.1 GHz. The result 
of such Planck spectroscopy is shown in Fig. 2(c). Clearly 
the quantum crossover temperature Tcr shifts to higher 
values with increasing frequency. Due to the finite un- 
certainty in Toff, we derive an effective crossover temper- 
ature Tcr + STcr [cf . red crosses connected by a red line in 
Fig. 2(c)] which again slightly deviates from the theoret- 
ically expected value Tcr = ^/2fcB [cf. dashed white line 
in Fig. 2(c)]. The magnitude STd- quantifies the mea- 
surement fidelity J" = 1 — \STcr\/Tci of our setup for 
vacuum fluctuations. Notably, for the entire frequency 
range > 95%. In summary, our Planck spectroscopy 
experiments not only provide clear evidence for the Bose- 
Einstein statistics of photons emitted by a conductor in 
the few photon limit, but also directly demonstrate the 



frequency dependence of the quantum crossover temper- 
ature. 

We next turn to the analysis of the mw beam splitters. 
Figures 3(a)-3(d) show the entire correlation matrix. The 
off-diagonal elements are cross-correlation functions mea- 
sured choosing (pi^o in order to obtain a maximum pos- 
itive result. This guarantees that the signals associated 
with the two detection channels are skewed in phase and 
no unwanted de-correlation is introduced during the ex- 
periments. Since the signal contributions of the thermal 
noise sources and the amplifier noise are independent, all 
cross-correlations vanish, e.g. (X1X2) ~ {xi){X2) = 0. 
Then, for = (3^ — 1/8 the covariance i?i2(0) = 
(^12(0) = cr^2 is obtained to 
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with the power cogain G12 = \[G\\fG2- We note that 
the temperatures Ta and Tq of the load resistors at port 
A and C are identical only in the ideal case, resulting 
in i?i2(0) = for the HR and WPD. However, in real 
experiments there may be small temperature differences 
ST = Ta — Tq. Figures 3(e)-3(h) show the covariance ma- 
trix as a function of T for w/27r = 5.0 GHz and using a 
WPD. The diagonal matrix elements iiii(O) and i?22(0) 
represent variance measurements and are analogous to 
the results of Fig. 2(b). The off-diagonal elements, in- 
stead, represent covariance measurements. It is evident 
that both the offset signal at 20 mK and the signal span 
between 20 and 350 mK for the covariance is reduced by 
about two orders of magnitude as compared to the vari- 
ance. This suggests that there is a cancellation of both 
the amplifier noise and the signal when measuring the 
covariance. The former is due to the fact that the ampli- 
fier noises are uncorrelated. The latter is expected from 
Eq. (2). In order to prove this conjecture, we use Eq. (2) 
to fit the experimental data. Doing so we use the cogain 
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FIG. 3: Full correlation function and covariance matrices measured at aj/27r — 5.0 GHz using a WPD (Gi ~ 90.5 dB, G2 — 
91.3dB, and G12 — 90.9dB. (a)-(d) Ruir) and Rijir) measured a,t T — 30 mK. (e)-(li) Temperature dependence of Rii{0) and 
Rij{0). The symbols represent experimental data, the lines fits to the data. 



which is determined by fitting the variance data. Fur- 
thermore, we set Ta — T, where T is the temperature 
measured by the thermometer, and use ST as the free 
fitting parameter. In this way, we obtain the red curves 
in Figs. 3(f) and 3(g), which are in excellent agreement 
with the data. We obtain 6T values of a few mK. From 
the covariance measurements we can get information on 
the characteristics of the WPD. Since Eq. (2) explicitly 
assumes the existence of four ports, the perfect fit of the 
experimental data provides clear evidence that the WPD 
effectively behaves as a four-port device. In the quan- 
tum limit, this fourth port adds vacuum noise to any 
given input signal. 

In order to confirm our findings on the WPD, we have 
measured the T dependence of the variance and covari- 
ance also for a HR (cf. Fig. 4) , which is a straightforward 
four-port beam splitter. The covariance data of Fig. 4(b) 
are in excellent agreement with the fitted curve obtained 
from Eq. (2). This clearly demonstrates that, both the 
HR and WPD are characterized by the same fundamental 
quantum-mechanical behavior. 

In conclusion, we have applied a correlation function 
analysis of the field quadratures to characterize black 
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FIG. 4: Temperature dependence of the (a) variance Rii{0) 
and (b) covariance Rij{0) measured at tj/27r — 5.85 GHz using 
a HR. The different power scale compared to Fig. 3(e)-(h) 
results from a different amplifier configuration. 



body radiation and the quantum properties of mw beam 
splitters. Our Planck spectroscopy experiments show 
that the mean thermal photon number emitted by a load 
resistor is following Bose-Einstein statistics and that the 
quantum crossover temperature shifts with frequency as 
Tcr = huj/2kB, providing an indirect measure of mw vac- 
uum fluctuations with high fidelity. Moreover, we have 
shown that, at the quantum level, vacuum fluctuations 
represent the minimum fundamental noise added by both 
three- and four-port microwave beam splitters. 
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